Coincidence detection of inharmonic pulses in a nonlinear crystal 
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Two trains of light pulses at periods that are equally shifted from the harmonics of a missing 
fundamental, are combined in a nonlinear crystal. As a result of a noncollinear phase matched 
second order nonlinear generation, a new train of pulses is obtained. When the temporal width of 
the input pulses is large, the frequency of the resulting pulse train follows the observations from 
classical experiments on the perception of virtual pitch by the brain. On the other hand, when 
the width of the input pulses is small, the generated pulse train exhibits much lower frequencies, 
analogous to those observed in the motor neural system. Our experimental set up allows us to 
explore, systematically and continuously, the transition between these two regimes, while at the 
same time demonstrate that coincidence detection in quadratic nonlinear systems has functionalities 
similar to those observed in the nervous system. 
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PACS numbers: 

One of the most remarkable features of nonlinear sys- 
tems is their ability to process complex input signals. A 
classical example is the perception of virtual pitch by the 
brain. In that context, it is well known [lj that a miss- 
ing fundamental tone can be perceived upon exposure to 
only some of its harmonics. Psycophysical experiments 
by Schouten et al. Q showed that when a set of con- 
secutive harmonics were equally shifted in frequency, it 
was not the frequency difference (i.e. the original funda- 
mental) that was perceived. Instead, the perceived pitch 
varied linearly with the frequency shift. Specifically, in 
the presence of input tones of frequencies 



fi = (fc + i -!)/„+ A/, 
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where /o is the missing fundamental frequency, k > 1 
is an integer and A/ is the frequency detuning which 
makes the input frequencies to be inharmonic. Under 
these conditions, the perceived pitch was seen to be given 
by a frequency that matches the following expression: 



fr = f + 
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Recently, Chialvo et al. Q proposed a simple and elegant 
mechanism that accounts for this response, involving a 
linear superposition of the input harmonics and a nonlin- 
ear noisy detection of the frequency (0 via a threshold. 
The mechanism, subsequently named ghost resonance, 
has been experimentally verified in lasers 0, 0, and in 
electronic circuits H- 
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Experimental results using magnetoencephalographic 
measurements Q showed that the missing fundamental 
illusion also arises when the harmonic inputs are pre- 
sented binaurally, i.e. different harmonics are applied 
to the each of the two ears. The mechanism of Chialvo 
et al. was extended to that situation by modeling 
separate neuronal pathways that detected two different 
input harmonics |9j- That study showed that in the 
context of distributed inputs, the mechanism of ghost 
resonance heavily relies of the coincidence detection of 
synaptic pulse trains (transduced by the input neurons 
that receive the input harmonic signals) by an integrat- 
ing neuron. An experimental realization of this effect 
in a real neurophysiological setup has been recently per- 
formed |lfj. That experiment has shown that when the 
input trains are inharmonic (i.e. frequency shifted with 
respect to the original harmonics), the processing neu- 
ronal pool responds at frequencies much lower that those 
expected from expression @. Subsequent experiments 
with nonlinear electronic circuits indicate that the 
difference in the response is due to the small width of 
the pulses acting upon the integrating neuron, in con- 
trast with the larger width presumably associated with 
acoustic neuronal pathways. 

In this paper, we use an utterly different experimental 
setup to address the question of the influence of the input 
pulse width on the ghost resonance response. This allows 
us to: (i) perform systematic measurements for continu- 
ously varying pulse widths, with high controllability and 
reproducibility, (ii) ascertain the generality of the phe- 
nomenon reported, which is seen to arise in any system 
that operates via coincidence detection upon threshold- 
ing, and (iii) suggest a possible functional role of this 
phenomenon in nonlinear photonic devices, with poten- 
tial applications in all-optical signal processing. 

The experimental setup is shown in Fig. 2] A 76 MHz 



2 



train of 130 fs pulses, produced by a Ti-sapphire laser at 
a wavelength of 800 nm, was divided in two beams using 
a 50/50 beamsplitter. These two beams are recombined 
in a BBO crystal cut for noncollinear phase matching, 
which generates second-harmonic light at a wavelength 
of 400 nm when both beams are present simultaneously 
in the crystal. The coincidence of the short laser pulses is 
obtained using a movable translation stage, as shown in 
Fig-d Both beams are chopped using electro- mechanical 
shutters Si and S2 at frequencies /1 and /2, respectively. 
The computerized chopping mechanism allows to control 
the shutter frequencies and ensures a constant (stable) 
phase relation between them. 
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FIG. 1: Schematic representation of the experimental array. 
We split the laser output, thin line, in two beams with a 
beamsplitter, BS, and determine the frequencies of the pulses 
using one shutter for each beam, SI and S2, both of them 
controlled by a code written in Lab VIEW (National Instru- 
ments). The thick dashed line is for the second harmonic at 
the exit of the nonlinear crystal, BBO; M are mirrors, F are 
filters to stop the pump beam, CC are corner cubes and MTS 
indicates micrometric translation stage. 

The input signals consist on two pulse trains of dy- 
namical frequencies f\ and /2, represented by the two 
upper thin traces in Fig. [21 In fact each pulse is, in turn, 
composed by a package of ultrashort pulses with a repe- 
tition rate of / rcp = 76 MHz and and optical frequency 
/opt = 3.75 x 10 14 Hz, as explained above (see inset of 
Fig. I2J). The width and dynamical frequency of each 
package are controlled by the optical shutters. In what 
follows we will refer to these packages as pulses, since the 
fast dynamics within the package is not relevant for the 
purpose of our experiment. 

When two input pulses coincide inside the nonlinear 
crystal (see Fig. |2J), a train of second- harmonic pulses at 
an optical frequency 2/ opt is generated by the crystal, 
and detected by the photomultiplier (lower pulse train in 
Fig. I2J . Therefore, the nonlinear crystal acts as a coinci- 
dence detector, and replaces the threshold of detection of 
previous systems where ghost resonance has been studied 
000 El In the particular case of Fig. [3 we have set 
fi = 2 Hz and f^ = 3 Hz, which corresponds to n = 2, 
k = 2, f = 1 Hz and A/ = in Eqs. Q-©. Given 
the thresholdless nature of the quadratic nonlinear inter- 
action, in this particular case a pulse train at the ghost 



frequency (f r — 1 Hz [see Eq. Q]) is always generated 
and its detection is only limited by the sensibility of the 
overall detection system. 
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FIG. 2: Time series of the input pulsed signals (two up- 
per lines, shifted vertically for clarity) of frequencies 2 Hz 
and 3 Hz, and of the response signal of the nonlinear crystal 
(lower line) at the ghost frequency 1 Hz, as detected by the 
photomultiplier. The inset shows the fully resolved time se- 
ries of a second-harmonic pulse, where the intensity does not 
go to zero because of limited response time of the detection 
system. Amplitude of the input signals has been vertically 
shifted in order to ease comparison between time series. 

In order to investigate the response of the crystal to 
inharmonic inputs and check the validity of Eq. |J2j) , we 
set the input frequencies to fi = kfg + A/ and fi — 
(fc + l)/o + A/, varying A/ between and 1 Hz. Since we 
are concerned about the influence of the pulse width on 
the pulse coincidence, we start by setting the pulse width 
At p to a relatively large value, namely At p = 60 ms. 
Figure Ofa) shows the instantaneous response frequency 
f r (defined as the inverse of the time interval between 
output pulses) as a function of slowest input frequency 
f\. The response frequency is seen to follow well the 
relation predicted by Eq. for k — 2 and k = 3 (dashed 
lines in the figure) for almost the whole range of /1 . Lines 
of different k are observed because the input frequency 
fx can represent different harmonics of the fundamental 
frequency /o, by changing the value of k. 

A linear response at large frequencies is also observed, 
however, in plot^a). This high frequency response is 
a direct consequence of consecutive coincidences of these 
broad pulses, as can be seen from Fig. H The resulting 
frequency is the inverse of the silent period of the slower 
input signal (T\ = -i- — Ai p ), as can be deduced from this 
figure, and therefore, for this case, the response frequency 
f r should follow the curve 
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This expression is represented by the upper thin solid 
line in Fig. |3Ja), and exhibits good agreement with the 
experimental results. 
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FIG. 3: Inter-pulse instantaneous frequency for increasing 
values of A/. The horizontal axis corresponds to the slow 
input frequency /i given by /i = fc/o + A/. The shutters were 
here driven for trains of rectangular pulses of width 60 ms 
(a) and 5 ms (b). The dashed lines in panels (a) and (b) 
indicate the expected ghost frequencies given by Eq. (|5J , and 
the upper thin full line in panel (a) represents the relation 
given in Eq. 



10 



0.04 



.0.02 



J 



3 4 
time (s) 

FIG. 4: Time series of the input signals (upper panel, shifted 
vertically for clarity) and response signal of the nonlinear 
crystal (lower panel), in the inharmonic case, for frequen- 
cies /i = 3.4 Hz and /2 = 4.4 Hz. We can see two domi- 
nant frequencies in this response: the lowest frequency corre- 
sponds to the values predicted by Eq. and the higher one 
arises from consecutive coincidences of the broad pulses used 
in these sequence of experiments. Here the pulse width was 
set to At p = 60 ms. 



We now turn our attention to the opposite limit of very 
narrow pulses. To that end we apply pulses of width 
Atp — 5 ms (the narrowest we can get given the speed 
limitations of the shutters). Figure E^b) shows the in- 
stantaneous response frequency f r , again as a function 
of slowest input frequency, /i . We can see that the rela- 
tion given by Eq. J2J) holds only in the neighborhood of 
the harmonic case (/i = 2 Hz and 3 Hz, f r — 1 Hz). On 
the other hand, in most of the inharmonic region, the 
system responds with pulse trains of very low frequen- 
cies, in agreement with previous experimental results in 
electronic circuits and the motor neural system [lCj . 
These responses are grouped in families of lines following 
f r = with b being an integer. The one with slope 1 
was observed experimentally in Ref. and its origin 
was determined analytically. The closer we are to the 
limit of zero-width pulses, the more lines appear. These 



lines arise from the condition of coincidence, 
I Ti=m T 2 , 



(4) 



where I, m are integers and the Ti and T 2 are the 
input periods, respectively, T\ = = k f \_^j an d 
T% = T£ = (fc +1 )^ 0+ A/ ■ When the condition given by Eq. 
Q is fulfilled, both input trains coincide and a pulse 
is detected. When potential coincidences at frequency 
A/ 7^ are themselves missed, lower frequencies —jf are 
generated. This gives rise to different families of lines 
[quasi- pyramids in Fig. EJ(b)], depending on the value of 
k that relates /i and /o- 

The two situations depicted in plots (a) and (b) of 
Fig.[2|represent two opposite limits of the coincidence de- 
tection mechanism and show that, in absence of a noisy 
threshold detection, the region of validity of Eq. (J2J) in- 
creases with the pulse width. In order to systemati- 
cally investigate the transition between them, we now fix 
the frequency shift A/ and vary continuously the pulse 
width. The corresponding result, for two different val- 
ues of A/, is shown in Fig. For A/ = 0.1 Hz, i.e. 
/i = 2.1 Hz and / 2 = 3.1 Hz (plot a), the k — 2 line 
appears only for pulse widths larger than 15 ms, and co- 
exist with slow frequency responses, which increase when 
the pulse width increases [ll|. The line corresponding 
to the case k = 3 does not appear for this detuning. 
For A/ = 0.7 Hz, corresponding to /i = 2.7 Hz and 
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FIG. 5: Inter-pulse instantaneous frequency for increasing 
values of A/ as a function of the pulse width. The input 
frequencies are (a) /i = 2.1 Hz and /2 = 3.1 Hz, and (b) 
/i = 2.7 Hz and fa — 3.7 Hz, which correspond respectively 
to a detuning A/ = 0.1 Hz (a) and A/ = 0.7 Hz (b) with 
respect to the harmonic k — 2 frequencies 2.0 Hz and 3.0 Hz 
[see Eq. @]. 

f% = 3.7 Hz (plot b), the response at k — 3 appears at 
pulse widths larger than 20 ms, and the one at k = 2 for 
widths larger than 40 ms. For broad enough pulses, the 
system responds at even higher frequencies. 

If we focus in the inharmonic region (when A/ ^ 0), 
the previous results clearly show the transition from a 
low-frequency response regime to another regime com- 
pletely dominated by relation as the width of the 
pulses increases. Such a transition between regimes could 
underlie the differences observed between previous psy- 
chophysical experiments on auditory response and re- 
cent experiments on the motor neural system |lOj |. In 
the former, expression J2J held unambiguously; in the 
latter, on the other hand, that behavior gave way very 
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frequently to low frequency responses. Correspondingly, 
synaptic pulses are known to be wider in the auditory 
system than in the motor reflex system. From a tech- 
nological viewpoint, the present results also show that 
nonlinear optical crystals carry out nontrivial signal pro- 
cessing tasks that mimic those of more complicated sys- 
tems such as the brain. 
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